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ABSTRACT 

In this report, the heat transfer and isothermal flow friction 
characteristics for seven compact heat exchanger surfaces are presented 
and compared. Five of the surfaces were of a similar triangular fin 
configuration, differing in material and method of construction, and 
two of the surfaces, having been tested previously, were of a modified 
rectangular passage configuration. The modified rectangular passage 
surfaces showed better overall performance than the triangular fin sur- 
faces, and matrices constructed from perforated material showed better 
overall performance than similar surfaces constructed from non-perforat- 
ed material. For similar matrix configurations, small variations in 
manufacturing procedure had little effect upon heat transfer and flow 
friction performance. 

The heat transfer data was obtained by means of the single-blow 
transient testing technique. A cylic method of transient testing was 
also investigated, and found to be both reliable and useful in that the 
additional method of testing extends the Reynolds number range for the 


experimental apparatus. 
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NOMENCLATU RE 


English Letter Symbols 


A Matrix total heat transfer area sq ft 
A, Base area (plane surface area) sq ft 
A. Matrix minimum free flow area sq ft 
Ae. Matrix total frontal area sq ft 
A Open area of matrix material due to sq ft 
perforations 
A Matrix solid cross sectional area sq ft 
§ available for thermal conduction 
A. Matrix solid cross sectional conduction sq ft 
area considering effect of matrix material 
perforations 
Cc Specific heat (gas) at constant Btu/(lbm deg F) 
P pressure 
C. Matrix material specific heat Btu/(lbm deg F) 
C. Flow stream capacity rate (mc) Btu/(hr deg F) 
C. Matrix heat capacity (Wicd Btu/ (deg F) 
Dy Hydraulic diameter of internal ft 
passage 
E Friction power expended per unit of hp/sq ft 
surface heat transfer area 
F Frequency of sinusoidal temperature cycles/sec 
oscillation 


G Matrix flow stream mass velocity (m/Ac) lbm/(hr sq ft) 


g Proportionality factor in Newton's 32.2(1bm ft) Albé a) 
second law 


h Thermal convection surface heat transfer Btu/(hr sq ft deg F) 
coefficient; heat transfer power per 
unit area per degree temperature difference 


K Loss coefficient for sudden flow contrac- dimensionless 
tion at matrix entrance 


Loss coefficient for sudden flow 
expansion at matrix exit 


Unit thermal conductivity 

Matrix material thermal conductivity 
Total matrix flow length 

Mass flow rate 

Pressure 

Matrix porosity (A./A,) 

Heat transfer rate 

Universal gas constant (53.35 for air) 
Electrical resistance 

Hydraulic radius (A_L/A,) 
Temperature 

Flow velocity 

Matrix volume 


Matrix material volume (corrected 
for effects of perforations) 


Matrix mass 


Distance from matrix inlet in 
flow direction 


dimensionless 


Btu/(hr sq ft deg F/ft) 
Btu/(hr sq ft deg F/ft) 
re 

lbm/hr 

lbf/sq ft 
dimensionless 

Btu/hr 

(ft lbf£)/(1bm deg R) 
ohm 

ft 

deg F 

ft/sec 

out t 


cu ft 


lbm 


ft 


Greek Letter Symbols 


p 


2, ¥ oP 


ca 1 <3 


Area density; compactness (ratio of matrix total 
heat transfer area to core volume-A/V_) 


Difference or change (time, distance, temperature) 
Fin efficienty 

Time 

Fluid viscosity (dynamic) 


Open area ratio due to matrix material 
perforations 


Density 
Ratio of free flow to frontal area (A/A,) 


Denotes "function of" 


Angular frequency of temperature oscillation 


sq ft/cu ft 


Sec, hr. 


lbm/hr ft 


lbm/cu ft 


dimensionless 


rad/sec 


Subscripts 


atm 


ave 


Local atmosphere 

Average 

Fluid (gas, air) 

Initial, inlet, individual 

Equitevent 

Mean, matrix 

At orifice 

Solid (matrix material), static 
Standard (temperature and pressure) 
Wall (solid surface) 

Local conditions 

At inlet (upstream of heating elements) 
At matrix inlet (downstream of heating) 


At matrix exit 
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Dimensionless Groupings 


f 


Ntu 


Fanning friction factor (Ratio of wall shear stress to the 
fluid dynamic head). This factor is plotted versus R to 
illustrate the flow friction characteristics of a matrix 
surface. 

*/3 
Colburn Factor C Nse Nev ), a generalized heat transfer 


grouping. This factor is plotted versus N. to illustrate heat 
transfer characteristics of a matrix surface 


mc 
Flow parameter ( Ws Cs ie) 
Nusselt number ( hOw/s ), a heat transfer modulus 
Prandtl number CA), a fluid properties modulus 
Reynolds number ( Di Yu ), a flow modulus 
Stanton number ( s/s ) heat transfer modulus 
GCp » a 


Number of heat transfer units ( h p% ), a heat transfer 


parameter 
f 


Ratio of matrix outlet temperature amplitude to matrix inlet 
temperature amplitude during cyclic temperature variation 


Ree er ce conduction parameter for solid matrix 
material ( "5 EVENT Cp mm 


ci Bearer baly conduction parameter for perforated matrix 
<sAk/, 
material ( “5 Uiilenee ys 


Distance parameter (x/L) 
Ti hA 

me parameter ( Ws Cs ) used in Locke's analysis 
Time parameter ( a Ye ) used in Cyclic technique 


Frequency parameter Cari ) used in cyclic technique 
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1. INTRODUCTION 

A compact heat exchanger may be defined as a heat exchanger having 
a large heat transfer surface area per cubic foot of core volume. 
Arbitrarily, heat exchangers with more than 200 square feet of heat 
transfer area per cubic foot of core volume are considered compact. 
Applications of particular interest for compact heat exchangers, are gas 
turbine power plants in aircraft, automobiles, and small seagoing craft. 
In all these applications the designer is limited by weight, flow frontal 
area, core volume, and pressure drop across the heat exchanger, An effec- 
tive basis for presentation and comparison of heat transfer and flow 
friction data for various heat transfer surfaces and core configurations 
is by plotting the dimensionless groupings- Colburn's heat transfer fac- 
tor, "j", and the Fanning friction factor, "f'', versus Reynolds number for 
heat exchanger surfaces. 

The data presented in this report was obtained by means of the 
single-blow transient testing technique described in the two sections 
following, and in reference [24]. The objectives of this thesis were 
to obtain and compare flow friction and heat transfer data for compact 
heat exchanger matrices of different configuration, to correlate this 
data obtained with data obtained for identical surfaces by means of 
a steady state steam-to-air testing technique, and to develop a testing 
technique that will extend the range of transient testing to higher 


Reynolds numbers. 
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2a SUMMARY OF THEORY 
A. BACKGROUND 
Analytical development of the theory for the single-blow or 
transient heating/cooling problem was performed in the 1920's by Anzelius 
(1926), Nusselt (1927), Hausen (1927-1929), and Schumann (1929). 
There was some duplication of effort in that Nusselt did not know 
of Anzelius' publication, and Schumann was unaware of the work of 
Anzelius and Nusselt. Anzelius and Schumann used similar procedures. 
Anzelius represented the gas and solid temperatures as unsolved, re- 
latively simple integrals, where Schumann derived and solved explicitly 
two infinite series utilizing Bessel functions. Schumann's solution was 
for a liquid flowing through a porous solid, however, dimensional analy- 
sis may be used to make his analysis applicable to gasses [24]. Hausen's 
derivation was similar to Schumann's, and in addition Hausen also solved 
the equations for cyclic heating and cooling. Nusselt's derivation was 
more simplified than Hausen's or Schumann's, but it is worth noting that 
Nusselt was the only one who considered the problem as part of a regenera- 
tor theory [10]. A summary of the single blow theory as it applies to 
a gas flowing through a porous solid is presented below. 
B. THEORY 
A homogeneous, porous solid has a fluid flowing through 
it. Both fluid and solid are at the same constant, uniform temperature. 
After a step change in the entering fluid temperature, the temperature of 
the solid and the fluid are to be determined as functions of time and 


position (along the flow length of the solid). 
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Starting with an energy balance in the matrix: 






mMCp (te x) 


hb(te -ts)dx 


eh 


lla al Ma 
EB 
“Kshs S38 Dolid 


b<——dx =o 







ts, att 
~ksAs ; ots.) 


Assumptions: 
(1) The fluid properties are temperature independent. 
(2) The fluid flow is steady. 
(3) The porous solid is homogeneous. 
(4) The thermal conductivities of the solid and the fluid are 
infinite in the direction perpendicular to the fluid flow. 
(5) The thermal conductivity of the fluid is zero in the flow 
direction. 
Boundary conditions: 
(1) The matrix is initially at uniform temperature. 
(2) At time zero, there is a step change in the entering fluid 
temperature. 


(3) The matrix boundaries are adiabatic. 
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Heat rates from the above energy diagram are: 


(1) Heat absorbed by the solid: PsAsCs oF dy 
(2) Heat transferred to the solid by convection: hb (t¢ -~ts)dx 
(3) Heat transferred from the fluid: MCp Eos 
(4) Heat transferred in the solid by conduction: - KsAs 2 OAS 1, 


The net resulting ater. balances for fluid and solid are: 
(1) Fluid: mep S SEE dx +hb (tr - igidn= 0 = (2-1) 
(2) Solid: P5As i SES > KsAs —ESux+hbitr-t s)dx (2-2) 


Now introduce the following dimensionless parameters: 





(1) Time parameter 
hA We 
Ts = Ke 2 
Ws Cs vr TE a 


Where: 


h = unit conductance for convective heat transfer 


A = matrix heat transfer area (sq ft) 
Wels = matrix capacity (Btu/deg F) 
O = time (hours) 
We = mass of fluid contained within matrix (1b m) 
Mm = fluid mass flow rate (lbm/hr) 
X = distance from matrix inlet in flow direction (ft) 


L 


Introducing the fluid specific heat, and regrouping, (2-3) becomes: 


te Dh hA xX \Wrp 


Ws Cs om Cp L Ws Cs 
Wer Cp 
3Cs 


is much greater than the thermal capacity of the fluid contained within 


matrix flow length (ft) 


(2-4) 


The term is very small since the thermal capacity of the matrix 


it. Thus, the dimensionless time parameter can be reduced to: 
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hA 
T= fava 9 (2-5) 





(2) Position parameter: 
a. hh 





WnCp fey r (2-6) 
Since the Number of Transfer Units, Neu = ao : 
z= Neu + (2-7) 
(3) Conduction parameter: 
KA 
= Se (2-8) 
mM apt 


Where; 
Ke = matrix thermal conductivity (Btu/hr sq ft deg F/ft) 


A == solid matrix cross sectional area &vailable for thermal 
conduction (sq ft) 


Introducing the above dimensionless parameters into the heat balance 


equations for the fluid and solid, and rearranging yields: 


Fluid: be” ts-te ” (2-9) 
Solid: obs Re AWN os (2-10) 


If thermal conduction in the solid is assumed zero in the direction 


of flow, the above equations simplify: 


gtr = ts-te (2-11) 


& 2-12 
Be = tr-ts ( ) 


From equations (2-11) and (2-12), Schumann, utilizing Bessel func- 


tions, developed his temperature-time-distance relationships: 


= on 
te-ti allen e. ee } a zn d( To (2iNtz)] 
tg mah’ -(2+T) a.) A int ) 
oF mer Sara < rte cul (2-14) 


aa dLt2)}" 


Plots of these functions versus 2 appear in Figure 1. Using these 
curves to determine Ntu from experimental data would be a tedious task 
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involving infinite series computations and curve comparisons. However, 
in his maximum slope technique, Locke [15] differentiated the theoretical 
curves, plotting slope versus ©/Ntu for various values of Ntu. Locke's 


expression for slope is: 
tp, -ti 
(pect) Neu“ (NE 
LF, -#) tu : "7 ut 
-J, (2iVTNtuwe end 
J (TAieu) (tT ku 








The downstream fluid temperature is taken at x = L so that 2 be- 
comes Ntu and te = Cee: Locke also evaluated the maximum slope of the 
temperature curve for each Ntu, and plotted maximum slope versus Ntu, a 
single curve. Figures 2 and 3 show Locke's curves. 

Both Schumann and Locke assumed thermal conductivity of the solid, 
in the flow direction, to be zero. The effect of longitudinal thermal 
conductivity is small in the high Reynolds number range, but it is con- 
siderable and should be taken into account in the low Reynolds range. 
To account for longitudinal thermal conduction, equations (2-9) and 
(2-10) must replace equ@€tions (2-11) and (2-12). Thus, the maximum 
slope will depend upon both Ntu and \ . Howard [9], utilizing a 
digital computer and finite difference technique, tabulated and plotted 
Ntu versus maximum slope for various values of \ . Howard's curves 
appear in Figures 4 through 6. By means of a properly designed experi- 
mental apparatus which satisfies Howard's idealizations and boundary 
conditions, the maximum slope of the heating or cooling curve and the 
conduction parameter can be determined. Ntu can then be found using 
Howard's data shown in Figure 5. 

In the vicinity of Ntu = 2, it can be seen from Figure 6 that there 
can be considerable error in Ntu for a small error in maximum slope. 
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In addition, it is more difficult to accurately determine maximum slope 
in the high Reynolds number region, as the maximum slope occurs earlier 
in time as Ntu decreases, and the determination of the maximum slope 
becomes difficult, having large uncertainties due to the inherent time 
response of the instrumentation. 

Another theory for determining the heat transfer characteristics of 
a porous solid, based upon a sinusoidally varying fluid temperature at 
the matrix inlet, was developed by Bell and Katz [2]. Starting as be- 


fore with an energy balance in the matrix: 





hb (tr -ts)dx 


“bk 


Solid 


— 


(1) All cross sections of the solid normal to the flow direction 


Assumptions: 


are uniform. 
(2) The fluid velocity is uniform throughout the matrix. 
(3) The fluid properties are temperature independent. 
(4) Thermal conductivity of the solid and fluid is infinite in the 


direction normal to fluid flow. 
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(5) Thermal conductivity of the fluid and solid is zero in the 
direction parallel to the fluid flow. 

(6) The matrix boundaries are adiabatic. 

(7) The convective heat transfer coefficient remains constant. 


Heat rates from the above heat balance are: 


(1) Heat absorbed by the solid: Ps AsCs se 
(2) Heat transferred to the solid by convection: hb (te-ts)dx 
(3) Heat absorbed by the fluid: ricp(Qte Awe: 3 Je) 
(4) Heat transferred in the solid by conduction: re) 


The net resulting energy balances for fluid and solid are: 


(1) Fluid: MmCp (bE dx + Bete) = hb (tr-ts) dx (2-16) 

(2) Solid: OsA5Co 258 dx = hb (t:-ts)dx (227) 
Inlet and initial — *: are: 

(1) Inlet: tp (0,6) = ty + ASIN WE (2-18) 

(2) Sinseiale eS) ee (2-19) 

(3) Initial: Ap &,o)= te (x) (2-20) 


where A is the amplitude of the sinusoidally varying temperature. 


Now introduce the following dimensionless parameters; 





(1) Time : T = Ae (2-21) 
(2) Distance: $= sis (2-22) 
(3) Frequency: w= wel (2-23) 
(4) Flow: Me rm Cp (2-24) 
“Ws Cs Glo 
(5) Heat transfer: hA (2-25) 
Néu 2 TIES 
Where: 
(A = flow velocity (ft/sec) 
GO = time (sec) 
L = matrix flow length (ft) 
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X = distance from matrix inlet in flow direction (ft) 


O&),= angular frequency of temperature oscillation (rad/sec) 


Introducing the above dimensionless parameters into equations (2-16) through 


(2-20) and rearranging yields: 


ote, dtr ee : 
(1) aks ot >t + Néu (tr ts) =O (2-26) 
(2) 268 -M,Ntu (tr-ts )=0 (2-27) 
(3) SeCore) = =n + ASIhH TT (2-28) 
(4) Cs (0) = es. (2-29) 
(5)" ae Cac eo (F) (2-30) 
This system has the periodic solution: 
-AG 
te (g,t)=ty+ Ae sih(wt-Gf) (2-31) 


-X5 
i aie CWC fm Néu sin (w -8§)- cos (wt -84)| (2-32) 


Where: Nétu 
~ bt (M Ntu)* 
MNtut 


< + W 
8 | 4 (iNtu)® 


A simple quantity to determine experimentally is R,, the ratio 





of the amplitude of outlet and inlet gas temperatures. Straightforward 


calculation shows: : - Niu 
Ae I+ (MNta)* 


A plot of Ntu versus Ra for various values of M is shown in Figure 7. 
The same experimental apparatus can be used for this cyclic operation 
as is used for the maximum slope method. The only addition needed is 
a method to sinusoidally vary the temperature of the air entering the 
matrix. 


Thus, heat transfer data can be determined for a wide range of 
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Reynolds numbers utilizing two testing methods on a single experimental 
apparatus; the single-blow technique producing heat transfer data for 


low Reynolds numbers, and the cyclic technique producing data for high 


Reynolds numbers. 


fay 


3% EXPERIMENTAL TECHNIQUES 
A. Single Blow Technique 


For Reynold's numbers up to approximately 1,000, the results 
of Howard's work [9] were used in determining the heat transfer data, 
Ntu and j, appearing in this report. In order to utilize Howard's 
theoretical data, the experimental apparatus had to conform to the ideali- 
zations stated in reference [9]. 

Idealizations: 

(1) The fluid flowing through the matrix remains steady and 
uniform in velocity and temperature at any cross section. The matrix 
thermal conductivity is infinite in the direction normal to fluid flow, 
and finite in the direction parallel to fluid flow. Thus the problem is 
one-dimensional in space. 

(2) The thermal capacity of the matrix is large compared with 
the fluid contained within it. This restricts the fluid to a gas, and 
means that there will be no time dependent terms in the equations for 
the fluid. 

(3) The thermal properties of the fluid and matrix are cons- 
tant and uniform. 

(4) The convective heat transfer coefficient is some suitable 
average and remains constant. 

(5) At time zero, the change in fluid temperature is a step 
change, with the matrix and its entrained fluid initially at uniform 
and equal temperatures. 

The test apparatus utilized for this report satisfies the above 


idealizations as described below: 
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(1) A square cross section contracting cone and screen type 
flow straighteners provided the necessary uniform velocity profile. In 
addition, the matrix test section was located as far upstream as possible 
to take advantage of this uniform velocity profile. The flow rate 
measuring apparatus was located downstream. The heater configuration, 
Figure 10, insured a uniform temperature profile due to even distribu- 
tion of heat over the flow cross section. The thermal conductivity in 
the direction parallel to fluid flow is taken care of by the use of the 
conduction parameter Y. 

(2) The use of a gas as the fluid flowing satisfies idealiza- 
tion (2). 

(3,4) A small temperature variation will insure uniform solid 
and fluid physical properties. At 80°F, a temperature variation of + 
10°F will result in an air viscosity variation of approximately + 1.5%. 
The changes in specific heats and thermal conductivities are negligible 
for a + 10°F change in temperature. 

(5) Various methods have been devised to produce a step change 
in matrix inlet temperature. The electric heater arrangement used in 
this report appears to produce the closest approximation to the desired 
step change with a minimum effect on flow properties [24]. 

A schematic diagram of the experimental apparatus is shown in 
Figure 8. Photographs of various components are shown in Figures 9 
through 14. The system is comprised of an entrance section with a contract- 
ing cone and a screen flow straightener, a heat section consisting of 28 
resistance type wire heaters made of nichrome wire in grid fashion, a 
matrix test section, an orifice metering section, a flow control valve 
and prime mover, and pressure and temperature measuring and recording 
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systems. Each of these components are further described in Appendix 
A. The experimental procedure for the single blow technique is as 
follows: 

Air is drawn through the system by means of a 30 horsepower Spencer 
Turboblower. The heater section is designed to heat the air to approxi- 
mately twenty degrees above ambient temperature. The mass flow rate is 
measured by means of an ASME orifice meter with removable orifice plates. 

Flow friction data is obtained by measuring the isothermal pressure 
drop across the matrix and the mass flow rate through the orifice meter. 

Heat transfer data is obtained as follows: the heaters are energi- 
zed to heat the air to approximately twenty degrees above ambient tempera- 
ture. When the matrix has been heated to a uniform temperature, the 
heaters are quickly de-energized and the time-temperature history of the 
air leaving the matrix is recorded on a Minneapolis Honeywell "Brown" 
strip chart recorder. 


The following data is recorded for each test run: 


(1) Atmospheric pressure Patm 
(2) Static pressure upstream of matrix | P 
(3) Pressure drop across the matrix AP. 
(4) Static pressure upstream of orifice Po 
(5) Pressure drop across the orifice AP. 
(6) Temperature upstream of orifice t. 
(7) Orifice diameter d) 


(8) Time-temperature cooling curve 


(9) Recorder chart speed 
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B. Cyclic Technique: 

To obtain heat transfer data by means of the cyclic tempera- 
ture variation method, the same testing apparatus can be used, the only 
modification being in the power input to the heaters. The idealizations 
set forth by reference [2] are the same as those imposed by reference [9] 
and shown on the previous pages with two exceptions: 

(a) The thermal conductivity of the solid is zero in the 
direction parallel to the fluid flow. 

(b) The temperature of the matrix inlet varies sinusoidal- 
ly. To satisfy idealization (a), this testing method is restricted to 
the high Reynolds number ranges ( Ne > 500), where Ntu is small (Ntu < 5), 
and the effect of longitudinal thermal conduction is small and can be 
neglected. (See Figure 15 and Chapter 5). 

To satisfy idealization (b), sinusoidally varying power is ap- 
plied to the heaters by means of a 1.5 KW Amplidyne receiving field excita- 
tion from a variable frequency function generator. This system is des- 
cribed further in Appendix B. 

To obtain heat transfer data, the temperatures of the matrix inlet 
and exit are recorded on the strip chart recorder. The amplitudes of 
the cyclic variation of these two temperature curves are then measured 
and Ntu is determined from the ratio of the amplitudes and the dimension- 
less parameter, M. Figure 15 shows a sample chart trace for a cyclic 
test. 
The same data are recorded as for the maximum slope technique. 
Data reduction relationships for both of these transient testing tech- 


niques are described in Appendix B. 
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4. DESCRIPTION OF MATRICES 

The surfaces investigated in this report are of two major varieties 
modified rectangular passage and triangular fin. 

The modified rectangular passage surfaces, previously examined by 
Bannon [1], and Piersall [22], are constructed of formed sheets separat- 
ed by flat splitter plates as shown below. The matrix approximates one 


made up of rectangular passages. The rectangular passages have an aspect 


ratio, Oo , of approximately 7: 





nickel and the other from .002 inch thick solid nickel. The perforations 
are of type 160/40 TV and are described further in Appendix C. | 
The five triangular fin matrices, provided by Solar, A Division of 
International Harvester, San Diega, were of varied material and method 
of construction. Geometric and physical data for these matrices appear 
in Figures 16 through 22. Two were made from nickel; one had fins made 
from .005 inch thick solid material, the other had fins made from type 
80/20 T perforated nickel material .0047 inches thick. The other three 
triangular fin matrices were made from .005 inch thick type 430 stain- 
less steel, and had no perforations. Of these stainless steel matrices, 
one had fins with very sharp bends, (all others had approximately sinusoi- 


dal corrugations), one had maximum fillet braze joints where fins and 
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splitters touched, (all others were not bonded by any. means), and the 
last stainless steel surface was similar in shape and construction to 
its solid nickel finned counterpart. All of the tgiangular fin matrices 
had a fin height of .1 inches and 10 fins per inch. All of the matrices 
measured 3-1/16 inch high X 3-1/16 inch wide. The modified rectangular 
passage matrices were 2 inches in flow length, while the triangular 
finned ones were 3 inches long. Geometric and physical data for all 


matrices appear in Figures 16 through 22. 
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5. PRESENTATION OF RESULTS 





Heat transfer and flow friction data for all surfaces tested are 
shown in TABLES I through VII and Figures 23 through 35. In Figures 
23 through 31 the Colburn heat transfer modulus, j, and Fanning fric- 
tion factor, f, are plotted versus Reynolds number. The Reynolds number 
was based on hydraulic diameter. In Figure 32, all surfaces are compared 
by a plot of j/f versus Reynolds number, and Figure 33 is a plot, for 
all surfaces, of heat transfer power versus flow friction power on a 
unit area basis. 

Figure 34 is a plot of Colburn j versus Reynolds number for the brazed 
stainless steel triangular fin surface. The j values for this plot were 
obtained by means of both the single-blow technique, and the cyclic tempera- 
ture variation technique described in Chapters 2 and 3. For Np < 1000 
the single-blow technique was employed, and for NR7 500 the cyclic tech- 
nique was used. This was the only surface tested by means of the cycli- 
cal temperature variation technique. 

Figure 15 shows a plot of Colburn j versus Reynolds number for the 
perforated nickel rectangular passage matrix. This plot shows the ef- 
fects of varying the longitudinal conduction parameter, X . One curve 
results from assuming N = 0, the second curve results from the 
calculated in reference [22] and described in Appendix C, and the third 
curve results from an electrical mn Logics calculation of XA which 
is also described in Appendix C. 

Figure 35 is a plot of Colburn j versus Reynolds number for the 
perforated nickel triangular fin matrice. the data for this curve was 
obtained by means of both the single blow technique and a steady state 
steam to air technique described in reference [29]. 
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In Figure 36, experimental and theoretical data are compared for 


the brazed triangular fin matrix. 
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6. EXPERIMENTAL UNCERTAINTIES 
The experimental uncertainties for the heat transfer and friction 
data appearing for this report were determined by the method described in 
reference [13] and shown in Appendix D. These uncertainties arise from 
possible errors in: 
(1) Physical constants - These physical constants were 
obtained from references [6], [7], and [8], and their uncertainties are 


estimated to be as follows; 


k, + -5% 
Cammmet 5% 
Come tes’% 
Ner + 2.0% 
P + 1.07% 
(2) Geometric measurements - These errors are due to con- 


struction inaccuracies and errors in linear measurement and are esti- 
mated to be as follows: 
Ac ,Am,A,As + 1-0% 
ul + .5% 
W, negligible 
(3) Instrumentation - Assuming adequate manufacturer's cali- 
bration a manometers and thermocouple wire, instrumentation errors will ' 
be assumed to be due to instrument fluctuations (manometers), and size 
of scale divisions (manometers and potentiometer). 
Pressures: 


Po +1.0% 


APyn + 1.0% 


P, + 1.0% 
AP, + 1.0% 


The error for the manometers depends upon the range of the pressure 
being read. In the range 1 to 3 inches of water, the error is .17% to 


.5% in the range 3 to 5 inches of water, the error is 1.0% to 1.7%. 


Temperatures: . 
fe tas deg Fe 


Maximum Slope: 
The measurement of maximum slope involves the drawing of 

a tangent and the linear measurement of three quantities xX , y , and 
m , as described in Appendix A. The percent error in maximum slope 
is estimated as + 2.0%. Entering gerd (7) with this value at Neu = 3.0 
and = 0, obtain Ntu error as + 7.2%. At Ntu = 25 and N= .05, Ntu 
error is + 10.0% for a 2.0% error in maximum slope. 

A summary of extreme values of uncertainty intervals calculated by 


the method shown in Appendix D is given below: 


m + 1.0% 
Ne + 2.3% 
J + 10% 


is 4%, 


late 
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7. DISCUSSION OF RESULTS 

Heat transfer and flow friction results appear in Figures 23 through 
36, with comparisons of all surfaces in Figures 31 through 33. 

In general, the triangular fin surfaces all had similar j and f 
curves; the nickel matrices having higher j and j/f values than the stain- 
less steel matrices. The high j and j/f for the nickel surfaces are due 
chiefly to the high thermal conductivity of nickel - 36.0 compared to. 
12.8 for 430 stainless steel (at 80°F). The see and j/f characteristics 
for the étatagsvee fin surfaces are in agreement with the theory presented 
in reference [16] for laminar flow through idealized cross sections (see 
Figure 36). 

The three stainless steel cores were quite similar in cross section, 
however, one might expect to find small differences in their heat trans- 
fer and flow friction characteristics due to slight differences in geo- 
metry. The rounded corners produced by the sine wave corrugations in 
two of the surfaces should reduce the effective surface area for heat 
transfer and wall friction as mentioned in reference [16]. The experi- 
mental results seem to agree in part with the above - the sharp cornered 
triangular passage geometry yielded higher a factors than the 
sinusoidal fin geometry, with the brazed sinusoidal fin core producing 
intermediate friction factors. The heat transfer characSeristics, how- 
ever, showed the j curve for the brazed core slightly higher than the j 
curves for the unbrazed ones. There are two reasons for this: First of 
all, the heat transfer area calculations for the brazed core took into 
account the loss of metal surface area due to brazing. In the unbrazed 
sinusoidal fin core, a portion of surface area similar to the "lost area" 
was relatively ineffective for heat transfer, but not accounted for in 
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calculating j. Secondly, the fin shape and matrix construction for the 
brazed core actually produced a geometry that is a cross between a tri- 
angular passage and a rectangular passage, with a rectangular passage 
having higher values of j and j/f than a triangular passage as shown in 
reference [16]. | 

It should be noted that these three stainless steel matrices yielded 
heat transfer and flow friction results that differed only slightly, in- 
aicating that small deviations in manufacturing procedure may have little 
effect upon heat transfer and flow friction performance. 

The rectangular passage cores had heat transfer and friction charact- 
eristics similar to the triangular fin cores; however, their heat trans- 
fer curves (j versus N,) were less steep, the perforated material pro- 
ducing a markedly higher curve than the solid material. The j/f flow 


area "goodness" factors for these surfaces increased as N. increased, 


R 
while j/f for the triangular fin surfaces decreased with increasing 
Reynolds number (See Figure 32). On the Dorp versus Eom plot, Figure 
33, rectangular passage surfaces showed definite superiority over the 
triangular fin surfaces, with the perforated material producing the best 
results. 

Figure 15 demonstrates the influence of conduction parameter upon 
the heat transfer modulus, j. Curves such as those shown dotted can 
result from an improper assumption ( AX = 0 in low Reynolds range), ap- 
proximations for ) (Appendix C), or an inaccurate value for Ko: 

Referring to Figure 7, it can be seen that there are areas where a 
small error in Amplitude Ratio, Ras would result in a large correspond- 
ing error in Ntu. However, the frequency of temperature oscillation, 


Je, can be chosen so as to produce a flow parameter, M, (M = me, NEC; 2) 
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oe ' that will result in Ntu being determined with minimum error for 
a given error in Rye For instance, at Ntu = 1.5, G),should be chosen so 
that M is less than .35, since values of M near .7 in this case result 
in a high uncertainty of Ntu. This flexibility is not inherent in the 
single-blow technique (where Ntu = 2). 

Figure 35 shows heat transfer data for the brazed triangular fin 
matrix. Both the single-blow and cyclic techniques were employed, thus 
producing data up to a Reynolds number of 1850 with an area of overlap 
between Reynoldd numbers 380 and 560. There seems to be fairly good 
correlation between the two testing methods with the exception of the 
data points at Reynolds numbers 450 and 560. These points were the re- 
sult of employing the single-blow technique in the area of high Ne and 
low Ntu (1.65 and 1.44), where Ntu is subject to error as explained in 
Chapter 2. The data from the cyclic technique was also subject to some 
error due chiefly to the cyclic power input to the heaters. As can be 
seen on the sample chart trace in Figure 37, there was some distortion 
in the temperature waveforms making them slightly different than true 
sinusoids. This distortion is due to inaccurate signal reproduction in 
the amplidyne. There was also a problem in instrumentation - the tempera- 
ture waveforms were recorded separately by means of the same instrument 
instead of simultaneously. Bell and Katz [2] utilized a dual channel 
instrument recording both temperatures simultaneously, reversing the in- 
put connections midway through each run so as to record each temperature 
on both channels, thus eliminating the effects of unequal channel sensi- 


tivities. The thermocouple configuration described in Appendix B does 


not allow for simultaneous temperature measurement and recording. 
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In Figure 35 data from the single-blow technique is plotted with 
data from a steady state stream to air technique. In both tests the 
identical type 80/20T perforated nickel triangular fin surface was in- 
vestigated. The steam to air test was conducted on a cross flow heat 
exchanger with the surface to be tested on the air side. The core con- 
figuration for the single-blow test is described in Chapter 4 and Figure 
17. Theoctitus in reference [28] utilized a dual testing technique similar 
to this in determining fin efficiency, Ue , and total surface temperature 
effectiveness te , for various surfaces. Theoclitus found for simple 
finned surfaces, that the theories for fin efficiency and total surface 
temperature effectiveness are reliable since he obtained a close corre- 
lation of data from the two methods of testing. For the steady state 
technique, theoretical values of Ne and 1. are used in determining the 
heat transfer characteristics, Ntu and j, while ": and I, have values 
of unity when the transient technique is employed, due to lack of a 
temperature gradient in the fin in the direction perpendicular to fluid 
flow (See Chapter 3). 

Figure 35 indicates a fairly good correlation of data for the two 
testing methods. 

Conclusions: 

(1) For matrices of the same geometry and material, small devia- 
tions in manufacturing procedure have little effect upon overall heat 
transfer and friction performance. 

(2) The single-blow and cyclic transient testing techniques 
utilized for this report are reliable in that the experimental results 


obtained agree with results obtained by means of steady state testing 
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and with theory. 
(3) The single-blow and cyclic techniques supplement each other 


in producing heat transfer data over a wide range of Reynolds numbers. 
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8. RECOMMENDATIONS FOR FURTHER WORK 
Refinement of the cyclic method of testing is needed both in theory 
and in testing apparatus. The theory could be extended to include the 
effect of longitudinal thermal conduction, this involving the use of 
another non-dimensional parameter in place of M, and a theoretical solu- 
tion using methods similar to Howard's [9] or Moreland's [20]. 
There are a number of areas for improvement in the test apparatus: 
(1) A modification is needed if higher flow rates are to be 
employed. Either a larger pipe diameter at the orifice 
meter, or another means of measuring flow rate is needed 
to increase the flow rates employed. 
(2) A more reliable power source is needed for the heaters so 
that a truer sinusoidal temperature variation results. 
(3) The temperature measuring system could be modified so 
that the matrix inlet and exit temperatures can be recorded 
simultaneously. 
It is also recommended that more surfaces be made available for 
testing by both transient and steady state methods so that data is avail- 
able for a wider range of Reynolds numbers, and the reliability of the 


two methods may be compared. 
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Figure 15. 
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APPENDIX A 
DATA REDUCTION RELATIONSHIPS 
GEOMETRY 


Porosity = Matrix free flow area 
Matrix frontal area 


Ac 
P = ee A-l 
Atr ot 
Hydraulic diameter = 4x Free flow area 
Wetted perimeter 


4AcL (A-2) 


Dy = mee 


Ap 


Area density = Heat transfer surface area 
Matrix volume 


A _ A 
8 ia Vin = Aey L (A-3) 


MASS RATE OF FLUID FLOW 
The mass flow rate, m, is calculated in accordance with references 
[21] and [23]. 


mM = 359k ee ae A NAP. ¥ (Ibm/hr) | (A-4) 


Where: 


k ra the flow coefficient including velocity of approach 
C = orifice coefficient of discharge from [21] 
{S = ratio of orifice diameter to inside pipe diameter 
O,= orifice diameter 
= thermal expansion factor 


= expansion factor 


= absolute static pressure at the orifice (1bf/sq ft) 


= 53.35 ft 1lbF, the universal gas constant for air 
lbm°R 
absolute temperature (°R) at the orifice 


Fo 
¥ = Ee the specific weight of the fluid flowing 
Pp 
R 
5 
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AP, = Pressure drop across the orifice in inches of water 


Making substitutions for kK and Y , (A-4) becomes: 


a 
Ter RT (A-5) 


In accordance with [23] figures (38) and (40b): 


Fa 1.0 


Y¥ 21.0 


Also: 
P = (PF ~ Fe ) .4912 X 144 (lbf/sq ft) 
am 13.6 : 
Where: woe is measured in inches of mercury 
Pe is measured in inches of water 
and: 


t= Co + 459.7 
where t, is measured in degrees Fahrenheit 
Introducing the values and expressions above into (A-5) produces 


the working formula; 


2 
4 aa ~ A-6 
we Ji - 8° om 


REYNOLD'S NUMBER 





By definition: 
- Dy? 
R= IM 
Substituting thes for Dy » and > for © i -: 
lo Cc 


(A-7) 


IN 
>: 


it 


| 


Ne 


:~ 
> 
= 
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FANNING FRICTION FACTOR 


UL LL LIAL LL LAA LLL LLM 

: 
a a 

C7227 ZAZA TTI TI TALIA 


TT PAL CLI Cee Poo Pte eee 


The following equation describes the flow through the matrix {12}: 


AP A 5 ENTRANCE ACCELERATION FRICTION Exi?t 
parte ouuM —~ pe Va Ab Vv _ b* VW 
P S| (ke pP) + 2(% 1) ee P'-K<e a | 


; Ac 
h = t 
where p = porosity ( her? 
It is assumed that UnpxVU, and Vx U2 since the associated 
pressure changes are very small compared with the total pressure. It 
is also assumed that the perfect gas low applies ( Us a ). Making 


the necessary substitutions and solving for ff produces for the iso- 


thermal case: 


Fs [2G OE - (RSP) (Eos) -(PB)(E-z)(rd]te ae 


RP 2 


Since the first term within the brackets above is much larger than 
the other terms, the following approximation can be made: 
P, + Pa ne 
——— = Py ~ I xPR 


Thus equation (A-8) becomes: 


£ 


2 Qe Gn a - (Ke+ Ke) + - (i+ e*) | = (A-9) 
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ke and en » the entrance and exit loss coefficients are functions of 
porosity, matrix configuration, and Reynold"s number, and are obtained 


in accordance with reference [11]. 


MAXIMUM SLOPE 
The maximum slope of the generalized cooling curve as previously 


mentioned, is a function of Ntu and conduction parameter, dM: 


‘ot 














d (TAnu) i —_ (A-10) 
Nee MAX 
where: 
_hA HA 
sc WsCs G and Neu = in op 
Therefore: 
TA = mice a8 
Ntu= To- 6 
. _ < 
and: d (TAiu) = a - 26 = Cr dO (A- 11) 
Cs 
alsoz te +: _ tr,- be, +] 
GF, -ti we -#} rapt” 
and: ae 
, EF. -£; es? ae te,-t; d (tr, -t) (A-!2) 
combining (A-11) and (A-12) yields: 
d (a | 
oeerse iy ee ae (tr, - ~ (A-13) 
d (Ayn) Fo tr,-6; | 
MAX 
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Now consider the following sketch of a generalized cooling curve: 





|: Gr, tn) i ¥ 
dO ie = 


doé=— 


CHART SPEED 


Y 


Mm 


d kr, -te,) 


Er, - ti) 


Combining the above produces the maximum slope: 


af cra) 
fia . Ss pee % ecHner SPEED 
d (en) Cr m 


MAX (A-14) 


Ntu can now be determined from conduction parameter » and maximum 
slope. 
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COLBURN  j 


Colburn's heat transfer parameter, j, is defined: 





) = 2/3 
J Nse Nor (A-15) 
= hh a 
Gep Nev “8 
= hA Ag 2/ 
mcp A Nev : 
thus: 
. 2/ 
j = Neu oe Ner ° (A-16) 
From equations (A-1) and (A-3): Ae A 
P= ahd coo 
Arr £ Aer L 
/3 1 PP 
= Né — — A-17 
From equation (A-17) it can be seen: 
id 
Ve R 
HEAT TRANSFER POWER 
The heat transfer power per unit area per degree temperature 
difference is [12]: 
Cp # : 
= J Ne “ 
h Nev Dn (A-18) 


For dry air at standard conditions equation (A-18) becomes: 


sto =-O2195 aa Btu 


. (he 2¢,F& degF) (4-19) 
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FLOW FRICTION POWER 


The flow friction power per unit area is defined as [12]: 
3 3 
| Vas { 3 
29c @C> \ Dy 7 Ne 
For dry air at standard conditions equation (A-20) becomes: 


= 
Ep = M10” (& )° F Nes) 


ney 
000 $¢. ft (4-21) 


For comparison purposes the surface geometries were reduced to a 


common hydraulic diameter. 
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APPENDIX B 
DESCRIPTION OF EQUIPMENT 
PRIME MOVER: 
Air as the working fluid is drawn through the system by means of 
a 30 HP, multistage, Spencer Turbo-Compressor. This compressor is de- 


signed for a flow rate of 550 cubic feet per minute. 


FLOW METER: 

Flow rates are measured by means of an ASME standard orifice section 
designed for convenient changing of orifice plates. Thin plate concen- 
tric orifices with throat diameters ranging from .308 to 2.31 inches are 
employed allowing for a wide range of flow rates. Pressure taps are 
located d upstream and d/2 downstream. Inside pipe diameter, d, is 


3.08 inches. 


HEATER SYSTEM: 
The heater section is made up of 28 grid type nichrome wire heaters. 
Two heaters are wound to each bakelite frame and thus may be energized 


in pairs. A schematic wiring diagram for one heater frame is shown below: 


L 
HEATCI2Z SYSTCM | HEATER 1, 










ON - OFF SELECTOR a 
SWITCH Mis aL { 
SS >... .., 
YOLTAGE . ee 


SOURCE 
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When energized, all heaters operate in parallel, the number of 
heaters in use depending upon flow rate. 

For maximum slope operation the voltage source is 220 volts AC and 
can be varied by means of a General Radio Company ''Variac" Autotrans- 
former. 

For cyclic operation the voltage source is a 1.5 KW General Electric 
Amplidyne connected in series with a D.C. source. The resulting input 
to the heaters is a sine wave voltage with an amplitude of 50 volts and 
an average value of 50 volts. Frequencies of .1 to 1.0 cycles per second 


are employed with a much wider range available. 


MATRIX HOLDER AND TEST SECTION: 
This section is constructed of polyethylene machined to close toler- 
ances. The T, and T, thermocouple grids and the matrix upstream and 


downstream pressure taps are mounted in the frame. The T, thermocouple 
grid is mounted in the removable matrix holder. 4%'' sheet styrofoam in- 
serts provide insulation and a snug fit for the test matrix. The flow 


channel measures 3-3/16"' X 3-3/16" and matrices up to 3" in length may 


be accommodated in the test section. (See Figure 12) 


TEMPERATURE MEASUREMENT: 

The temperature at the orifice meter is measured by means of a 
copper-constantan thermocouple and read in millivolts with a Rubicon 
Portable Precision Potentiometer. 

All other temperatures are measured by means of iron-constantan 
thermocouples, each temperature being measured by a group of five 
thermocouples connected in series. The Tt thermocouples are bound to- 


gether and insulated from each other by teflon tape. This bundle is 
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mounted at the throat of the intake cone, upstream of the wire screen 


flow straightener. The T T,,and T, grids are connected to a multiple 


i 1 


barrel switch so that the emf generated by T, may oppose the emf generat- 


1 


ed by either To or T.- Thus the difference between ambient and matrix 
inlet or exit temperature can be read. These temperature differentials 
are recorded on a Minneapolis-Honeywell "Brown" Strip Chart Recorder with 
variable sensitivity and chart speed. For maximum slope operation, T, 


versus T is recorded. For cyclic temperature variation, Ty versus T> 


and Ty versus T are recorded. 


PRESSURE MEASUREMENT: 

The four pressure taps shown in Figure 8 are connected to Ellison 
draft gages by means of Imperial Company "Poly-flo" plastic tubing. The 
draft gages are a combination of inclined and vertical gages. They are 
shown in Figure 9. For small pressure differentials, Vernon Hill Company 
Type "C'' Micromanometers are used. Pressures of 1.25" to 3.0" of water 
are read by means of the inclined draft gages. Pressures below this range 
are measured by means of the micromanometers; Pressure above 3.0" are 


read on vertical draft gages. 
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APPENDIX C 
CONDUCTION PARAMETER FOR PERFORATED MATERIAL 
The porous material used in the perforated plate matrices investi- 
gated in this report was made by Perforated Products, Inc. by means of 
an electro-depositing process. 


For a solid plate or fin, the conduction parameter is defined as: 


_ KsAs 
d Le 


where A. is the solid cross-sectional area for conduction which is 
identical to the product of plate width and thickness, and L is the 
conduction path length which is identical to the matrix flow length. 

For a perforated material such as the type 160/40 TV shown in 
Figure 38(a), the cross-sectional area for conduction varies in the 
flow direction, and the conduction path length is greater than the flow 
length. 


The definition of conduction parameter is based upon a solid materi- 





are. ksAs5 
= LC  RyC 

where Ren = sz the matrix thermal conduction resistance. 
s s 


Using the above definition for a perforated material, an "equivalent" 
conduction parameter, dk ,» can be evaluated by determining the thermal 


resistance (R By defining a relative shape factor S', the conduc- 


th k’ 


tion parameter of a perforated material, XK ; can be related to the con- 


duction parameter of the same material as if it were unperforated: 


rv ikKsA 
i= Sipe eS ircies 


Or: 


Rr, = 2 ( Ron) 
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By the analogy between conduction heat flow and electric current, note 


that for similar shapes, the ratio of electrical resistances: 


ee = 5S! 
Te (Ri) 


In Figure 38(b) are shown two shapes, each representing the shaded 
portion of Figure 38(a). The shape on the left represents an unperfor- 
ated material, and the shape on the right is an incremental model for the 
160/40 TV material. The relative shape factor is determined from the 
ratio of the electrical resistances of these two models. Since the 
material is of uniform thickness, two dimensional models can be construct- 
ed from electrically conducting paper (such as Teledeltos paper) as shown 
in Figure 38{c), with low resistance silver paint applied as indicated. 
The electrical resistances are then measured with an ohm-meter. This 
analog method produced the following results for the type 160/40 TV per- 
forated nickel material: 


| RT, he = 8.5 = 
Sic RT 27 SIS) 


thus: he = ,3IS 

Whereas, Nie = .127 as determined in reference [22] using: 
= \ Be we 

Ng =) As Lig 

A.= > An; 


Lye ZL; 
and: ALS and Les are approximations indicated in Figures 38(a), and 


(b). 


where: 


The effects of these approximations upon Ntu and j are shown in 


Figure 15. 
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APPENDIX D 


EXPERIMENTAL UNCERTAINTY INTERVALS 


According to reference [13], the uncertainty interval for a 


result is: 
aR .,\*. /aR a OR 1% 
= kc wi) + (35 Wa) Stee Vi, Wn) | (D-1) 


Wy» Wos++-W are the uncertainty intervals of the factors making 


up the result. 


R is the expression for the result 


V 


V peoeV are the factors making up the result. 


ae? 

The expression for Colburn j is: 
2 

Neu - Npy (s 


The uncertainty interval for Colburn j is: 





2 (a 
Wj si | @ Te Wuew) + + (2b Wace) f DA wp) + (Qe. 2/5 WN pr? “s) | 
4 


“a Cr Nee vs) i (ie New va) (Nis AcNpr's Wa) + (3 Heese ed] 


WiC Mee + ley +0 + MEY] 


d 


For Ntu = 3.0; A= 0 etc: 
Wi = 7.5% 
N) 


For Ntu - 25.0; \#= .05 ete: 


aie 19.27% 
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APPENDIX E 
DIGITAL COMPUTER PROGRAM FOR DATA REDUCTION 

Included herein is a computer program for a CDC 1604 digital com- 
puter. By means of this program all necessary calculations were made 
to produce the heat transfer and flow friction data appearing in the 
tables and graphs of this report. 

A data reduction program was written by Bannon [1]. This program 
was then further sophisticated by Piersall [22], producing conduction 
parameter » » maximum slope, and friction factor as well as other flow 
information from geometric and experimental data inputs. 

The present program will now produce all necessary heat transfer and 
flow friction data in table form from a minimum of input. When data from 
the single-blow technique is processed, a curve-fitting interpolation 
subroutine is utilized to determine Ntu from maximum slope and conduction 
parameter. The input for this subroutine is Howard's [9] maximum slope 
data (Table VIII). For the cyclic technique, Ntu is determined from 
| equation (2-33) by means of an iteration process. 

The program will process data obtained by means of either the 
single-blow or cyclic technique and data from more than one matrix can be 


reduced in each computer run. 
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SUBROUTINE CYCLE (RB, EM, P) 


COMMON RB, EM, P 


P = 0.0 
62 RC = EXPF(-P/(1.0+(EM*P)**2) ) 


IF (RC-RB)64, 64,63 


63 P = P + .002 
60 TO 62 

64 RETURN 
END 
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